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ABSTRACT 

We present new mock catalogues for two of the largest stage-IV next-generation 
surveys in the optical and infrared: LSST and Euclid, based on an N-body 
simulation+semi-analytical cone with a posterior modification with PhotReal. This 
technique modifies the original photometry by using an empirical library of spectral 
templates to make it more realistic. The reliability of the catalogues is confirmed by 
comparing the obtained color-magnitude relation, the luminosity and mass function 
and the angular correlation function with those of real data. 

Consistent comparisons between the expected photometric redshifts for different 
surveys are also provided. Very deep near infrared surveys such as Euclid will provide 
very good performance (Az/(1 + z) ~ 0.025 — 0.053) down to H ~ 24 AB mag and 
up to z ~ 3 depending on the optical observations available from the ground whereas 
extremely deep optical surveys such as LSST will obtain an overall lower photometric 
redshift resolution (Az/(l + z) ~ 0.045) down to i ~ 27.5 AB mag, being considerably 
improved (Az/(1 + z) ~ 0.035) if we restrict the sample down to i~24 AB mag. Those 
numbers can be substantially upgraded by selecting a subsample of galaxies with the 
best quality photometric redshifts. We finally discuss the impact that these surveys 
will have for the community in terms of photometric redshift legacy. 

This is the first of a series of papers where we set a framework for comparability 
between mock catalogues and observations with a particular focus on cluster surveys. 
The Euclid and LSST mocks are made publicly available. 

Key words: catalogues - surveys - galaxies: luminosity function, mass function 
galaxies: photometry - galaxies: evolution - cosmology: observations 


1 INTRODUCTION 


In the next decade, several large-area cosmological surveys 

are expected to be completed in the optical: the Javalam- 
bre Physics of the Accelerated Universe Astrophysica l Sur¬ 
vey (J-PAS, iBem'tez et al.ll2009al : iBenitez et alJ 20141 ). ex¬ 
pected to start in 2015; the Dark Energy Spectroscopic 
Instrument (DESI, iLevi et ahl l2013i ) survey, expecting to 
start in 2 01 8; the Large Synoptic Sky Telescope (LSST , 
llvezic et al.l 120081 : ILSST Science Collaboration et aD 1 20091 ) 
survey, also starting in 201 8; t he near infrared: the Eu¬ 
clid survey (ILaureiis et al.l l201ll). which will start in 2020 


and the Wide-Field Infrared Survey Telescope (WFIRST0), 
which will start in 2020 as well, among others. One of the 
main goals of such surveys is the measurement of Bary- 
onic Acoustic Oscillations (BAOs) up to redshift 1.3, there¬ 
fore achieving a combined figure of merit (FoM) superior 
of 100 0, making them to be clas sified as Stage IV experi¬ 
ments ll Albrecht et al.l[2006l . 120091 ). Different strategies have 
been designed to achieve such goals: J-PAS will be sam¬ 
pling the optical spectrum with 54 narrow-bands, conse¬ 
quently obtaining a superb photometric redshift resolution 
(Az < 0.003(1 + z)) for all galaxies down to flux limit of 
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i ~ 23.5 mag (iBenitez et al.l[2014h : DESI will target 20- 
30 million galaxies and quasars spectroscopically over 14000 
degrees square; LSST will image a large portion of the sky 
with 6 broad-band optical bands down to very faint mag¬ 
nitudes (r ~ 27.5), consequently obtaining very deep qual¬ 
ity data; Euclid will obtain very deep infrared space data 
(H ~ 24 mag), providing a extremely deep and good reso¬ 
lution high-redshift data sample and WFIRST will obtain 
very deep photometry (J ~ 27 mag for 2400 square degrees 
and H ~ 29 mag for 3 square degrees) in six broad bands 
in the infrared, sampling also the high-redshift universe. All 
these surveys will be covering thousands of degrees of the 
northern and southern sky. This, together with the collec¬ 
tion of external datasets (X-rays and Sunyaev-Zel’dovich ef¬ 
fect maps, spectroscopic samples, etc) will result into an 
enormous quantity of available data, several order of mag¬ 
nitudes superior to the present data in terms of combined 
area, depth and photometric redshift resolution. 

Up to date, different studies have been performed to es- 
timat e the photometric redshift performance of different sur- 
veys Ovaizu ct al. 2008; Bancrii ct al. 2008; Benitez ct al. 
2009a; Hildebrandt et afl 120101 : jB^agamba_et_alj_ 2012|; 


Dahlen et al.l 20131 : Gorecki et ahll2014l : Molino et al.l 12014 ) 


with a variety of different techniques. Two approaches have 
been considered in these works. The former uses empirical 
spectroscopic samples to calibrate the photometric sample 
while the latter extracts their results from simulations. The 
first methodology, while using real data, provides usually op¬ 
timistic results since the spectroscopic samples are usually 
biased towards brightest samples. Also, these methods pro¬ 
vide usually an average global accuracy for the whole popu¬ 
lation of galaxies at different redshift and magnitude ranges, 
making these results difficult to be compared consistently. 
The second methodology allows a treatment of spectroscop¬ 
ically complete mock samples, hence allowing to provide ex¬ 
pected results as a function of magnitude or redshift. 

However, it is well known that semi-analytic galaxy 
formation models are not fully representing the observa¬ 
tional universe. Some of the inconsistencies between those 
models and the observations are related to the inconsis- 


tency with the stellar mass function 

(Fontanot et al. 20091; 

Hansson et al. 20121; Skelton et al. 2012; Weinmann et al. 

20121: Henriciues et al.| 2013, EoiH 

Mitchell et al. 20131), 

luminosity function (Fontanot et al. 

2009: Cirasuolo et al. 

20101; Somerville et al. 2012|; 

Skelton et al. 20121: 

Hansson et al. 2012: Henriaues et al. 20121; [Guo et al. 2013a!; 

Gonzalez-Perez et al. 2014: Henriaues et al. 201311. chemical 

abundances (Fontanot et al. 2009i: 

De Lucia & Borgani 


Skeltonet^L 201^_ [Somerville et alJ ~ 20121 : Hansson et al.l 
20121 : iHenriaues et al.l l2012fl . consequently producing a 


wrong tilt of the color-magnitude relation, ’plume’ effects of 
redder galaxies belonging to the cluster spread within 1-2 
bright magnitudes, absence of a smooth-transitory green 
valley and other related effects. These effects generally 
lead to an underestimation o f the photometric r edshif t 
uncertainties . As shown by lArnalte-Mur et all (l2014h : 
lAscaso et al.l (120151 1 . obtaining photometric redshifts di¬ 
rectly fr om moc k ph otome try mimicking the ALHAMBRA 
survey dMoles et ali 120081 4 provided a mean photometric 
redshift dispersion t hree t imes high er than the expected 
for the real data ((Molino et al.112014 1. After creating a new 


photometry set using a technique that uses a new empirical 
library of templates (see §2.2) to fit the original spectral 
distribution, the photometric redshift accuracy exactly 
matched that obtained from real data. 

In order to create a solid basis for comparison between 
different methodologies and datasets, we start with this pa¬ 
per a series of works under the title: Apples to apples, A 2 , 
which intend to provide a handy well-characterized mock 
catalogue with realistic well-behaved photometric-redshifts 
that mimic two of the next generation surveys: LSST and 
Euclid. A separate work (Ascaso et al. in prep) will also be 
devoted to analyze the multiple narrow-band J-PAS survey. 

In this work, w e have a pplied a similar pres cription as in 
lArnalte-Mur et all 120141) ; (Ascaso et al.l (1 201 5l ). The whole 
point of a applying this prescription is to produce a galaxy 
sample which resembles as much as possible the real universe 
in terms of galaxy properties, in particular their distribu¬ 
tion of luminosities, redshift, and spectral types. The cho¬ 
sen mock catalogue is ba sed on the pu b licly available light 
cone mock catalogues bv IMerson et al.l (I2013h . These cata¬ 
logs are extracted fr om an N-body simula tion from the Mil¬ 
lennium simulation, l(Springel et al .1(200,51) and se mi-analytic 
model of galaxy forma tion from GALFORM (I Cole et al.l 
120001 : iBower et al.ll2006l) . We have use d the Bay esian Photo¬ 
metric Redshift software 2.0 (BPZ2.0. iBenftez (l2000l L Ben¬ 
itez et al 2015, in prep), to obtain photometric redshifts. 
This technique, called PhotReal (Benitez et al. in prep), is 
chosen to provide really accurate and realistic results at com¬ 
paring with real data and avoid fake effects such as ’plume’ 
or unrealistic color-magnitude relations. 

At least two more works will continue these A 2 series. 
In the second one (Ascaso et al., in prep), we will pro¬ 
vide consistent clusters selection functions for all the sur- 
veys considered by using t he Bayesian Cluster Finder (BCF, 
lAscaso et al]|2012l . 12 014al ). a well behaved detector which is 
maximizing the advantages of many methods in the liter¬ 
ature to detect not only red-sequence cluster and groups. 
The third A 2 paper (Ascaso et al., in prep) will be devoted 
to forecast the cosmological parameters from the obtained 
cluster counts in the same consistent way. These mock cat¬ 
alogues will be publicly available and they will be useful 
for any consistent prediction between the considered sur¬ 
veys: LSST and Euclid in terms of galaxy evolution, galaxy 
clustering, large-scale structure or any other astrophysical 
purpose. 

The structure of the paper is as follows. In section 2, 
we introduce the mock catalogue used in this work and we 
give an overview of the technique developed to provide re¬ 
liable photometry and realistic photometric redshifts. Sec¬ 
tion 3, summarizes the main characteristics of the three 
next-generation surveys considered in terms of area, depth, 
number of bands and different technical properties. Section 
4 provides an extended comparison with different observa¬ 
tional dataset to check that their main properties such as 
stellar masses, magnitudes, colors and spectral types resem¬ 
ble those obtained from the PhotReal post-processed mock 
catalog. Section 5 describes the obtained results on photo¬ 
metric redshift resolutions and the fraction of outliers for 
each survey. Section 6 provides a way to improve the over¬ 
all photometric redshift performance by obtaining the ’best 
selection samples’ with an odds quality selection and finally, 
section 7 summarizes the results of this work and discusses 
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the possible applications of the overall set of data to future 
studies on providing cluster selection functions and cluster 
cosmology constraints as well as the impact of this science. 

The mock catalogue used in this paper has a fixed cos¬ 
mology of fiju = 0.25, 12 a = 0.75, f lb = 0.045 and h = 0.73, 
chosen to match the cosmological paramet ers estimated 
from the first year results from the WMAP (iSpergel et al.l 
120031 ) and which are not significantly different to t hose ob¬ 
tained by Planck (IPlanck Collaboration et al.ll2013l ). Hence, 
throughout the whole paper and when necessary, we will 
use the same cosmology. All the magnitudes in the paper 
are given in the AB system. 


2 MAIN SIMULATION 

2.1 Original light mock catalogue 

In this work, we based our s i mulat ions on the original 
mock catalogs by iMerson et al. | (j2013tPl The spatial infor¬ 
mation of this mock catalogue comes from the Dark Matter 
halo merger trees ext racted from the Millennium Simulation 
dSpringel et al.ll2005l ). a 2160 3 particle N-body simulation of 
the AC DM cosmology. This simulation traces the cold dark 
matter structures from redshift z=127 through the present 
day in a cubic volume of 500 h~ 1 Mpc 2 . The minimum reso¬ 
lution of this simulation is 20 particles, corresponding to a 
halo resolution of 1.72 x10 10 /i _1 Mq. 

These haloes have been populated with galaxies cre - 
ated with the variant of the GALFO RM dCole et al .1 120001 1 
semi-analytical model bv I Bower et al.l (1200611 . which models 
many physical processes governing the formation and evo¬ 
lution of the galaxies such as the star formation and merger 
history, the radiative cooling of the gas, the feedback as a 
result of supernovae active galactic nuclei, etc. It also makes 
predictions for many galaxy properties including luminos¬ 
ity as a function of wavelength. However, these predictions 
do not always resemble those of the observed galaxies: the 
colors of the faint galaxies tend to be too red compared to 
the real ones, there exist a ’plume’ of redder galaxies than 
the red sequence at redshift lower than 1.4, etc. We have 
overcome these problems by running PhotReal (see i )2.2l) . 
an algorithm that creates more realistic photometry by us¬ 
ing a set of empirical templates to fit the original spectral 
distribution. Finally, these galaxies have been associated a 
sub-halo according to its status of central or satellite galaxy 
at each snapshot. More inf ormation on the con struction of 
this mock can be found in iMerson et al.l (120131 ) and refer¬ 
ences herein. 

We have chosen to use the Euclid 500 deg 2 wide mock 
catalogue since it is the deepest of the surveys considered 
and hence, we can always make it shallower. As for the sur¬ 
veyed area, we chose to use a 500 square degrees mock even 
if it is a factor of 20 to 40 times smaller than the total 
area of the targeted surveys as a compromise between hav¬ 
ing a easy-to-transfer, relatively small in size catalogue and 
an area large enough to have significant enough statistical 
results for our scientific purposes. Indeed, since the main 
purpose of A 2 is to analyze photometric redshift and clus¬ 
ter selection performance, none of these results should be 


2 http://community.dur.ac.uk/a.i.merson/lightcones.html 


substantially affected after applying a correction which in¬ 
creases the errors following a Poissonian approach to the 
results to account for this effect. 

As for galaxy clusters and groups, the possible struc¬ 
tures that we might be missing due to the limitation of the 
area would be the rarest ones (M > 10 15 Mq ) with similar 
properties as those recently discovered such as, fo r instance 


the SPT-CLJ2344-4243 (The Phoenix, McDonald et al 


2012), ACT-CL J0102-4915 (El Gordo , iMenant eau et al 
2012fh SPT- CL J0546-5345 jBrodwin et al.ll2010l) . SPT-CL 
J2106-5844 dFolev et al.ll201lh at redshift > 1. These struc- 
tures, which have been reported to be very unfrequent 
djee et al.ll2009l : iHovle et aid 1 20121) . are usually the easiest 
to detect among other structures at the same redshift al¬ 
most with any technique (optical, X-rays, cosmic shear) so, 
we can ensure that if they exist, we will be detecting them. 
The other possible problem that might arise from using a 
smaller simulation would be related with those structures 
lying on the edges of the cone. However, when detecting 
galaxy clusters and groups in the next paper or these se¬ 
ries, we will use the BCF which, automatically corrects the 
probabilities by the missing area of the clu ster due to edge 
or masks issues dAscaso et alJl20l3 . l2014al ) . 


2.2 PhotReal: the photometry calibration 
technique 

In order to obtain more realistic and reliable photometry 
and photometric redshifts, we have a pplied PhotReal , a 
similar algorithm as in prev i ous works (Amalte-Mur et al .1 
I20l4 IZandivarez et akll2014l : lAscaso et ah 20151 ) to the ini¬ 
tial mock catalogue. We summarize here this technique de¬ 
veloped to fit the spectral types to make the photometry 
more realistic. A more extended explanation will be given 
in a forthcoming work (Benitez et al. 2015, in prep). 

This technique relies on a library of empirical templates 
developed for the well-known photometr ic redshift cod e, the 
Bayesian Photometric Redshift fBPZ2.0. lBenftej|2000l . Ben¬ 
itez in prep). This library consist of eight empirical tem¬ 
plates with a very low outlier rate and bias in high quality 
photometric catalogues. Hence, it contains a complete (up to 
1-2%), even if coarse-grained, representation of real galaxy 
colours for the galaxy populations sampled. 

The PhotReal procedure starts by running BPZ2.0 in 
’ONLY_TYPE yes’ mode (similarly to BPZ) to match the 
optical ugriz and infrared (if available) rest frame photom¬ 
etry of the initial mock catalogue to the best template from 
the empirical library. In other words, PhotReal forces the 
semi-analytical photometry to have a realistic template dis¬ 
tribution by choosing the closest BPZ2.0 template. Once the 
SED is thus chosen, we generated galaxy flux es th ro ugh t he 
survey set of filters in the sam e way as in the iBenftez et al.1 
(2009allbh: lBenitez et al.l d2014l) empirical mocks and add to 
them empirically calibrated photometric noise. This noise is 
a combination of the expected photometric noise from the 
observed relationship between magnitudes and errors in a 
similar characteristics survey filters (see §3), plus a system¬ 
atic noise which is approximately constant with magnitude 
and most likely unavoidable when measuring galaxy colours 
in multiband photometry. This systematic is empirically cal¬ 
ibrated to be 8% for bluer objects and 6% for red galaxies 
(Benitez, in preparation). 
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Afterwards we run BPZ2.0 on those mock catalogues to 
obtain the final photometric redshift estimation and spec¬ 
tral type. This procedure of course implicitly introduces an 
error: there may be real galaxies which are present in the 
catalogues, but not covered by the BPZ2.0 templates or not 
adequately present in the spectroscopic redshift catalogues 
used to measure the outlier rate. Since we are not aware of 
any such substantial population within the depths for the 
surveys considered, we estimate that the procedure followed 
here does not introduce a significant bias. 

The final PhotReal output includes photometry and 
photometric errors, spectral types, stellar masses, a single 
photometric redshift estimation and the redshift probability 
function P(z,T), condensing the probability as a function 
of redshift and spectral type, for all galaxies in the mock. 


3 CONSIDERED NEXT GENERATION 
SURVEYS 


Traditionally, the estimation of the accuracy of photometric 
redshifts has been given globally, without specifying their 
dependence wi th photometry or redsh i ft (although see re- 


cent efforts by Hildebrandt et al. 2012j; 

Dahlen et al.l 2013!: 

Raichoor et al. 20141: Molino et al. 12014 

); it has been esti- 


mated for different spectroscopic samples with very com¬ 
plicated selection functions to reproduce and different flux 
cuts have been performed. These issues carry a number of 
arduousness. First, it becomes complicated to have a clear 
idea about the best region of the parameters space to use if 
one wants to select the best quality photometric redshifts. 
Also, these assumptions make very hard the comparison be¬ 
tween the photometric redshift accuracy of different surveys. 
Furthermore, a clean sample obtained by selecting the best 
quality photometric redshift has been proved to have many 
benefits for many scientific purposes but the previous points 
produce a non equivalent treatment for different surveys. 

In the spirit of A 2 , we pretend to perform a systematic 
comparison between the photometric redshift estimation for 
different next-generation stage IV surveys using the same 
mock catalogue and the same photometric redshift tech¬ 
nique. In this section, we highlight the main characteristic of 
the two next-generation surveys we have considered in this 
work in terms of design, area coverage, wavelength range, 
photometry depth, complementary data and other technical 
issues. 


3.1 LSST 

The L arge Synoptic Survey Telescop <J^|( LSST , llvezic et al .1 
120081 : ILSST Science Collaboration et ai1l2009l l is a stage IV 
experiment starting in 2018. The main objectives of this sur¬ 
vey are also the constraints of the dark energy parameters, 
the weak lensing analysis and the study of transients. 

The survey is designed to be a large, wide-held ground- 
based survey imaging the whole visible southern sky from 
an 8.4m telescope Cerro Pachon (Chile) with six broad-band 
optical bands ugrizy down to r=27.5 mag after coadding 


10 years of operations. The depth of the survey will al¬ 
low to produce very detailed shear and weak lensing maps 
in a simil ar way that was done with the Deep Lens Sur¬ 
vey (DLS, Wittman et alJl 2002 j ; ISchmidt fe Thormanll2013l : 
lAscaso et al.l 2014al). T he exp ected depths can be found in 
Table 1 in llvezic et al] (12008 1. 


The photometric errors have been estimated accord¬ 
ing to the description given by the Survey Science Group: 
eLSSlQ, e.g. the expected error for a galaxy in a single ob¬ 
servation is 


4 4 I 6 

& sys ■ ®photom 

with a 3ys is the systematic photometric error ranging be¬ 
tween 0.005-0.01 mag for the LSST case, and a p hotom is the 
random photometric error, which can be written as 

Ophotom = (0-04 - 7 )* + ■yx 2 , with x = io 0 ' 4(m_m5) 

where ms is the 5 a limiting magnitude depth and 7 depends 
on the sky brightness and readout noise. Eventually, the final 
error obtained for N observations is a n = 


3.2 Euclid 

The Euclid surve\Fl (|Laureiis et al.l l201l] ') is an European 
space mission survey starting in 2020. Euclid will also be 
a Stage IV experiment and it is expected to provide very 
accurate measurements of the dark energy constrains from 
a number of different probes including BAOs, weak lensing, 
cluster counts and high-z supernova. 

The survey will have a ’Euclid Wide Survey’ covering 
15,000 square degrees of the sky in the near infrared (Y JH), 
down to ~24 mag in H band providing also near-infrared 
spectroscopy and two ’Euclid Deep Fields’, about 2 magni¬ 
tude deeper than the wide survey which will cover around 
20 square degrees each. The deep infrared imaging will map 
and systematically explore the mostly unknown-high red¬ 
shift Universe. 

Euclid space observations will be combined with other 
space and ground-based observations to obtain the source 
photometric redshifts and physical properties. Among the 
optical surveys that will be available form the ground, 
one is already available, the Sloan Digital Sky Sur¬ 
vey (SSD S), and two are planned, the LSST an d the 
DES (iThe Dark Energy Survey Collaboration! l2005ll sur¬ 
veys. Hence, we will consider two cases which will provide us 
a maximum and minimum limit in the photometric redshift 
and cluster survey predictions for the Euclid survey: 

• Euclid pessimistic case: where the optical data coun¬ 
terpart will come from the five broad band photometry 
gr izY from t he DES down to the depth stated in Table 1 
bv lMohr et al.l (l2012l l. The photometric errors for the DES 
band s have been e stimated from the mock catalogues by 
IChang et al.l (l2014f) which mimic the properties of the DES 
data. 

• Euclid optimistic case: where the optical counterpart 
will come from the combination of DES and LSST-like sur¬ 
veys. 


3 http://www.lsst.org/lsst/ 


4 http://ssg.astro.washington.edu/elsst/magsfilters.shtml 

5 http://www.euclid-ec.org/ 
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In both cases the photometric errors for the ’Euclid 
Wide survey’ have been calibrated from the Cosmic Assem¬ 
bly Near-inf rared Deep Extragalactic Legacy Survey (CAN- 
DEL^fl IGuo et all 1 201 3bll and norma lized to the Euclid 
depth (e.g. H~24 at 5a, Laureiis et al. B). 


4 PERFORMANCE OF THE MOCK 

GALAXIES 

4.1 Model properties of the galaxies 

With the double-fold objective of measuring the variation of 
the modeled properties of the galaxies used in this simula¬ 
tion and proving the non-dependence of the results shown 
here on t he photometric redshift code used, we have used 
LePhare dllbert et al.l I2006I . l2009l l to obtain an estimation 
of the stellar masses, ages and star formation rates (SFRs) 
of the galaxies in the original mock and the PhotReal post- 
processed catalogues. 

To do so, we have used 27 stellar population synthesis 
models generated using lBruzual fc Chariot! (12003 1 models at 
14 ag es th at ran ge fro m 0.1 to 13.0 Gyr. These BC03 models 
use a IChabriei] I2OO3I) initial mass function with metallici- 
ties of Z=0.004, 0.008, and 0.2 Zq and are characterized by 
exponentially declining star formation histori es (S FH s) with 
timesc ales of r = 0.1 — 30 Gyr. We adopt the lCalzetti et ahl 
ll2000l l reddening law with 0 ^ E(B — V) ^ 0.5. We first 
obtained the photometric redshift estimation and we later 
fixed it to provide model measurement of the galaxies. 

Then, we have compared on one hand, the LePhare stel¬ 
lar masses estimated from the mock initial photometry with 
the original stellar masses in the mock and on the other, 
the LePhare stellar masses estimated from PhotReal post- 
processed photometry with the stellar masses obtained with 
BPZ2.0. In Fig. [l] we display the difference between the 
LePhare estimation to the original photometry and the ini¬ 
tial mock catalogue masses (solid line); and the difference 
between the LePhare estimation to the PhotReal photom¬ 
etry and the stellar masses estimated with BPZ2.0 to this 
photometry (dashed line). The difference between the initial 
catalogue masses and LePhare estimations is a factor of two 
more biased towards negative values (< SM(LePhare)i — 
SM{Mers<yn) >= —0.028 ± 0.494) than the difference be¬ 
tween the BPZ2.0 masses and their estimation with LePhare 
(< SM{LePhare)f - SM{BPZ2.0) >= -0.018 ± 0.432). 
However, the dispersion obtained in the relation is similar. 

Looking into this issue more in detail, we have plotted in 
Fig. [2] the density map of the logarithmic ratio between the 
stellar mass given in initial mock catalogue and its fit with 
LePhare as a function of the initial stellar mass for different 
redshift bins. While the agreement of the stellar masses at 
lower redshifts is tight (0.114), the dispersion increases as 
we increase in redshift, finding that the masses estimated 
with LePhare are systematically smaller th an the original 
ones. A similar behavior has been observed in lMitchell et al.l 
(12013P when comparing the initial input masses with mod¬ 
eled masses and it has been associated to the presence of 
dust in the model. 


http: //candels.ucolick.org/ 



log(M'(LP)/M‘(Sim)) 

Figure 1 . Solid line: Distribution of the difference between the 
original stellar masses provided in the original mock catalogue 
and the stellar masses fitted with LePhare to the original mock 
photometry. Dotted line: Distribution of the difference between 
the stellar masses obtained from the PhotReal post-processed 
mock catalogue and the stellar masses fitted with LePhare to the 
PhotReal mock photometry. The vertical dashed line indicates 
the zero value for reference. 


Similarly, we have plotted in Fig. [3j the logarithmic ra¬ 
tio between the stellar mass obtained from the fit to the 
PhotReal photometry and the estimation with LePhare as a 
function of the PhotReal mass as a function of redshift. In 
this comparison, we see that in the lower redshift bin, the 
bias of the stellar mass estimation is smaller than the previ¬ 
ous comparison (-0.002). However, the dispersion is 3 times 
larger. As we increase in redshift, we n otice a bi modal distri¬ 
bution of the difference, explained in iMitchell et al] d2013l l 
for the use of different Initial Mass Function (IMF) and 
Stellar Population Synthesis (SPS) models. However, the 
dispersion do not increase, becoming comparable between 
0.5 < 2 < 1.0 and a factor of 1.5-2 smaller at z>1.0. 

As a consequence, the PhotReal post-processed cata¬ 
logue seems to be a factor of ~ 1.5-3 less sensible to the 
presence of dust in the model. In 44.51 we will investigate 
more in detail the stellar mass function. 

In Fig. H we show the differences between stellar ages 
and star forming rates (SFR) obtained from the initial and 
transformed mock catalogue for different redshift bins. The 
stellar ages of the considered galaxies are very similar, par¬ 
ticularly at high redshifts, where they fully agree. Overall, 
we find that more than ~ 60% of the sample has a difference 
in galaxy ages less than 0.5 dex. 

As for the SFR, we find that the differences are smaller 
in the highest redshift bin (1.5 < z < 3), obtaining < 
log(SFRi/SFRf) >= —0.281 ± 0.509. At lower redshift 
the dispersion increases a factor of 3-4 due to the presence 
of a secondary peak at all redshift bins with a systematic 
difference of log(SFR) < —3. These secondary peaks indi- 
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Figure 2. Density map of the logarithmic ratio between the stel¬ 
lar mass given in initial mock catalogue and its fit with LePhare as 
a function of the initial mock stellar mass for six different redshift 
bins. 
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Figure 3. Density map of the logarithmic ratio between the stel¬ 
lar mass obtained from the fit to the PhotReal photometry and 
the estimation with LePhare as a function of the PhotReal mass 
for six different redshift bins. 


cate that a substantial part of the galaxies in the PhotReal 
post-processed mock are better fitted by a spectrum with 
different star forming rates than originally defined, partic¬ 
ularly at lower redshifts. This is compatible with the shift 
to bluer colors shown in the previous section where we show 
how the PhotReal post-processed catalogue makes the color- 


Figure 4. Normalized distribution of the difference in logarithm 
between galaxy ages (left column) and star forming rates (right 
column) obtained with LePhare from the initial photometry and 
the PhotReal post-processed photometry. From top to bottom, 
the differences are shown as a function of increasing redshift: 0 < 
2 < 0.5, 0.5 < 2 < 1.0, 1.0 < 2 < 1.5 and 1.5 < 2 < 3.0, 
respectively. 


magnitude relation (CMR) more realistic showing a blue 
cloud as expected in observational data. 

While the objective of this test is not to obtain reli¬ 
ably model parameter of the galaxies in the mock but rather 
check the difference between their properties, we conclude 
that the main star formation history of the galaxies in the 
mock are not changing substantially, particularly at high 
redshift. 


4.2 Spectral type recovery 

The original mock catalogue bv lMerson et al.1 (120131 ) did not 
include the information regarding the spectral type of the 
galaxy or its spectrum due to disk space storage. We have in¬ 
stead performed a comparison of the spectral type obtained 
from the fit to the original photometry, ti, and the recovery 
spectral type obtained from the PhotReal photom etry, tf. 
The spectral types are similar to the ones used in iBenfted 
(l200Cll) ; the first time refers to a typical Elliptical, type 2 to 
a Lenticular, type 3-4 to Early and late spirals respectively, 
and type 5-8 to starburst. In Fig. 0 we show these type 
distributions for the Euclid-Pes (top left panel), Euclid-Opt 
(bottom left panel) and LSST (top right panel) surveys re¬ 
spectively. 

The difference between the original and recovered spec¬ 
tral type distribution is minimum for the Euclid-opt case 
(1.90%), as expected and it changes slightly for the Euclid- 
pes case (3.65%) and the LSST case (2.90%). While a higher 
departure of the final spectral distribution is observed for the 
LSST and Euclid-pes, we conclude that the spectral distri¬ 
bution is well-recovered for the three cases considered. 
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Figure 5. Spectral type distribution obtained from the original 
Merson photometry (£j, solid line) and from the PhotReal post- 
processed mock catalogue ( tf , red dashed line) for each of the 
considered surveys: Euclid-Pes (top left panel) and Euclid-Opt 
(bottom left panel) and LSST (top right panel). The bottom panel 
shows the difference between the respective distributions. The 
Euc-opt shows an excellent recovery of the spectral type, whereas 
the difference between the original and the recovered simulation 
for the Euclid-pes and the LSST becomes a factor of 1.5-2 larger. 

4.3 Color-magnitude relation 

In the original semi-analytical mock catalogue, the first ef¬ 
fect that we noticed is the fact that the CMR usually ex¬ 
hibited a lack of blue cloud that it is not observed with 
observation a l data (e.g. lAscaso et al . 20081: iMei et al., 200l 
lAscaso et al.ll2014bl) in agreement with lMerson et al.l d2015h . 
PhotReal overcomes this problem by re-computing the pho¬ 
tometry and therefore, shifting some galaxies from the red 
sequence to the blue cloud. Note that the magnitudes that 
we recover have already incorporated photometric error es¬ 
timated for the survey whereas the original mock catalogues 
do not. 

In order to illustrate this effect, we displayed in Fig. [6] 
the g — i colo r versus i magnitud e density plots resulting 
of stacking the iNelson et all (l200ll ) CMR in redshift bins of 
0.1. The CMR information has been extracted using DexteiQ 
and the original V — I versus I relation have been converted 
to g — i versus i magnitude by adding a color term. On top 
of these density maps, we overlapped, in the four top pan¬ 
els, the CMR contours obtained from stacking the original 
(four top panels) and PhotReal (four bottom panels) post- 
processed photometry using the LSST bands for the same 
redshift bins. In all cases, the contours limit are shown to 
5% of the total density CMRs. 

Similarly, we plotted the CMR density maps at redshift 
2 > 0.8 in Fig [3 The top left panel displays the r — z versus 


7 http://dexter.sourceforge.net/ 



Figure 8 . z — J color versus J magnitude diagram for a random 
case of a 2.17 X10 14 M cluster at z = 0.52 in the mock catalogue as 
seen by the Euclid-opt survey with its bands. The left panel refers 
to the original mock CMR whereas the right panel represents the 
recovered CMR after introducing the real effects of the catalogs. 
The blue squares show the ’pl ume ’ galaxies in the original CMR 
and their location in the PhotReal post-processed CMR. 

z magnitude at the redshift slice 0.8< z <0.95; the top right 
panel shows the i—z versus z magnitude for the redshift slice 
1.1 < z < 1.3 and finally, the bottom left panel refers to the 
z — J versus J magnitude for the redshift bin 1.55 <z <1.65. 
The first two CMRs are extracted from IMei et ahl (120091) , 
also with Dexter, whereas the last CM R correspo nds to t he 
galaxy clusters at redshift 1.62 by IPapovich et ahl (|2010l ). 
These colors are used since they straddle the 4000 Abreak 
in each redshift range. At the view of these two figures, we 
observe that PhotReal reproduces the colors of galaxies in 
the blue valley, in agreement with the observational data, 
whereas the original mock catalogue displays a lack of blue 
galaxies. 

A secondary effect related to colors, noticed in the orig¬ 
inal mock catalogues is the presence of a ’plume’ of redder 
galaxies than the CMR. While this effect has been noticed in 
the infrared colors at high-redshift (z^l.4), the effect contin¬ 
ues appearing at lower redshift in the optical bands. In Fig 
[H] we show an example of the z — J color-magnitude relation 
of a cluster at z ~ 0.52 and M ~ 2.17 x 10 14 Mq as seen by 
the Euclid-Opt survey in the original mock catalogue (left 
panel) and in the PhotReal post-processed one (right panel). 
The blue squares indicate the galaxies considered to be in 
the ’plume’ in the original mock catalogue and their position 
in the PhotReal post-processed mock catalogue. We observe 
that the presence of this systematic ’plume’ disappears from 
the color-magnitude diagrams after the correction but they 
are still red galaxies, located in the main red sequence of the 
cluster. These galaxies have not changed substantially their 
colors and, in addition, they have also diluted their photom¬ 
etry within the photometric errors of the simulation. This 
result ensures us that, while we are making the mock cata¬ 
logue more realistic, we are not changing substantially the 
properties of their galaxies. 

4.4 Luminosity Function 

An additional check that we have performed is the compar¬ 
ison between the infrared and optical luminosity function 
(LF) at all redshift ranges. In Fig. [9] we plot the observa- 
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Figure 6. Stacked g — i color versus i magnitude density maps from Nelson et al. 2001 for different redshift bins (0.30< 2 : <0.40 top left 
panel (585 galaxies), 0.40< 2 <0.50 top right panel (713 galaxies), 0.50< 2 <0.60 bottom left panel (430 galaxies), 0.65< 2 <0.80 bottom 
right panel (296 galaxies) respectively for the four top and bottom panels). The color scale is logarithmic and redder colors represents 
higher densities while bluer colors refer to less dens regions. We overplottted the density contours of the stacked CMR obtained from 
the mock photometry using the DES and LSST bands for the same redshift slices. The four top panels display the contours obtained 
from the initial photometry, whereas the four bottom panels refer the contours obtained from the post-processed PhotReal photometry. 
Phot Real reproduces in an accurate way the galaxy color distribution in both the red sequence and the blue valley within the redshift 
range 0.30< 2 <0.80. 
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Figure 7. Stacked r — z color versus 2 : magnitude density map for the redshift bin 0.8< 2 : <0.95 (332 galaxies); i — z color versus 2 : 
magnitude density map for the redshift bin 1.1< 2 : <1.3 (224 galaxies), both redshift bins from Mei et al. 2009 and z — J color versus 
J magnitude for the redshift bin 1.55< 2 : <1.65 (266 galaxies) from Papovich et al. 2010 (left top, right top and left bottom panels 
for the three top and bottom panels respectively). The color scale is the same as in Fig. [6] We overplotted the density contours of the 
stacked CMR obtained from the mock photometry using the DES, LSST and Euclid bands for the same redshift slices. The three top 
panels display the contours obtained from the initial photometry, whereas the three bottom panels refer the contours obtained from the 
post-processed PhotReal photometry. PhotReal reproduces in an accurate way the galaxy color distribution in both the red sequence and 
the blue valley within the redshift range 0.80< 2 : <1.65. 






















10 Ascaso, Mei & Benitez 



Figure 9. H-band rest-frame luminosity function for three differ¬ 
ent redshift bins: 1.5 < z < 2.0 (top panel), 2.0 < z < 2.5 (middle 
panel) and 2.5 < £ < 3.0 (bottom panel). The red diamonds refer 
to the magnitude counts of the initial mock catalogue and the blue 
triangles refer to the same counts for the Phot Real post-processed 
mock catalogue for the Euclid-opt survey. The solid line and the 
gray area are the observational luminosity function obtained by 
Stefanon & Marchesini, 2013. 

tional H-band rest-frame LF for thr ee different redshift bins 
at z > 1.5 obtained from the fits bv I Stefanon &; Marchesinil 
(I2013I) from the Multiwavelength Survey by Yale-Chile 
(MUSYC, iGawiser et al.l l2006h . the Fai nt Infrared Ex- 
tragalactic Survey (FIRES, iFranx et al.l [2000) and the 
GOODS Chandra Deep Field-South Survey (FIREWORKS, 
IWuvts et al.ll2008h . We overplot the luminosity counts ob¬ 
tained for the Euclid-opt simulation in the original mock 
catalogue and in the PhotReal post-processed one. We find 
a very good agreement between both counts from the simu¬ 
lation and the observational LF at least down to Mh < —21 
where the flux limit of the simulation is being reached. This 
becomes a further indication that the photometry in the H- 
band at least it is not significantly altered. 

We have also checked the concordance of the LF at lower 
redshift. In Fig. [TUI we show the i LF for three d ifferent red- 
shift bins as measured by iMontero-Dorta fe Pradal (120091 ) 
for the SDSS DR6 survey at low redshift (0< z <0.1) and 
iGabasch et al.l (l2006h from the FORS Deep Field (FDF) sur¬ 
vey for 0.45< z <1.3. As before, we overplot the luminosity 
counts of the initial and PhotReal post-processed LSST sim¬ 
ulation for all the three redshift bins. The counts obtained 
for both mock catalogues are comparable and reproduce well 
observational luminosity function at any redshift. 

4.5 Stellar mass function 

We have looked into the agreement of the observational stel¬ 
lar mass (SM) function with the one obtained from our simu¬ 
lations. In Fig. 1111 we display the observational SM obtained 



Figure 10. Rest-frame i luminosity function for three different 
redshift bins: 0. < z < 0.1 (top panel), 0.45 < z < 0.85 (model 
panel) and 0.85 < z < 1.3 (bottom panel). The red diamonds 
refer to the magnitude counts of the initial mock catalogue and 
the blue triangles refer to the same counts for the PhotReal post- 
processed mock catalogue for the LSST survey. The solid line and 
the gray area are the observational luminosity function obtained 
by Montero-Dorta & Prada (2009) for the low redshift bin and 
Gabasch et al. 2006 with the Case 3 evolution for the other higher 
redshift bins. 


bv lMuzzin et al.l (I 2 OI 3 I) from the COSMOS/UltraVista field 
within 0 < z < 3 for six different redshift bins. As before, 
we overplotted the SM counts for the original mock pho¬ 
tometry and the PhotReal post-processed photometry. We 
find a remarkable good agreement between the observational 
SM function and both estimations in the mock catalogues, 
giving also support to the fact that we do not substantially 
alter the initial mock catalogue. 

We have also compared the SM function for different 
galaxy types. I n Fi gs. [TUI and [TU I we show observational SM 
functions bv lMuzzin et al.l (120131) for quiescent galaxies and 
star forming galaxies respectively, we also display in these 
plots the SM counts for these types of galaxies for the orig¬ 
inal and PhotReal post-processed catalogue. To do this, we 
have defined quiescent galaxies (star-forming galaxies) as 
those that the spectral classification obtained from the fit 
to the original L, or PhotReal post-processed photometry, 
tf has a spectral type <3 (^3). 

In general, we find a very good agreement between the 
observational SM function for the star-forming galaxies and 
the initial and PhotReal post-processed SM distribution for 
any redshift bin. There exists also a good agreement between 
the observational SM function for the quiescent galaxies and 
both mock SM, at least for redshift higher than 1. Interest¬ 
ingly, the high mass end of the mass function obtained with 
PhotReal reproduces better the observational stellar mass 
function, compensating for the fact that the original mock 
catalogue shows lack of quiescent galaxies in this mass range. 
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Figure 11. Stellar mass function for six different redshift bins: 
0.2 < z < 0.5 (top left panel), 0.5 < z < 1.0 (top right panel), 
1.0 < z < 1.5 (middle left panel), 1.5 < z < 2.0 (middle right 
panel), 2.0 < z < 2.5 (bottom left panel) and 2.5 < z < 3.0 
(bottom right panel). The red diamonds refer to the stellar mass 
counts of the initial mock catalogue and the blue triangles refer to 
the same counts for the PhotReal post-processed mock catalogue 
for the Euclid-opt survey. The solid line and the gray area are the 
observational SM function obtained by Muzzin et al. 2013. 

4.6 Angular correlation function 

As a final test for our mock catalogues, we have also com¬ 
pared the obtained angular correlation functions with obser¬ 
vational data. For this purpose we used the public software 
athen4E Since the basis of our cosmological simulation is 
the Millennium simulation and we have not perturbed the 
galaxy positions, the overall correlation function should re¬ 
flect reality. However, our aim is to test if the angular cor¬ 
relation function of passive and star-forming galaxies has 
unrealistically changed when applying PhotReal. 

In Fig [H we display the overall angular correlation 
function for the initial and the PhotReal post-processed 
mock limited to redshift 0.2. This is compared t o the angular 
corre lation functions obtained from the SDSS (IWang et al.l 
120131) split for different magnitude bins. For the sake of clar¬ 
ity, we have not split our sample in magnitude bins. How¬ 
ever, we see that our estimation is well located within the 
different observational angular correlation functions. 

Moreover, we have split the angular correlation func¬ 
tions for quiescent and star forming galaxies with the 
same criteria for classifying the initial and PhotReal post- 
processed photometry as in the previous section. In Figs. 
El and HU we display the angular correlation function for 
the original and Euclid-Opt PhotReal mock catalogue for 
the quiescent and star-forming galaxies respectively for six 
different redshift bins. We observe that the angular corre¬ 

8 http://www.cosmostat.org/athena.html 



Figure 12. Stellar mass function for quiescent galaxies for six 
different redshift bins: 0.2 < z < 0.5 (top left panel), 0.5 < z < 
1.0 (top right panel), 1.0 < z < 1.5 (middle left panel), 1.5 < 
z < 2.0 (middle right panel), 2.0 < z < 2.5 (bottom left panel) 
and 2.5 < z < 3.0 (bottom right panel). The red diamonds refer 
to stellar mass counts of those galaxies classified as early type 
(it < 3) in the original mock catalogue fit. The blue triangles 
refer to the stellar mass counts of those galaxies classified as early 
type (it < 3) in the PhotReal post-processed mock catalogue for 
the Euclid-opt survey. The solid line and the gray area are the 
observational quiescent SM function obtained by Muzzin et al. 
2013. 

lation functions for the both types of galaxies are almost 
indistinguible for the original and PhotReal post-processed 
catalogue, confirming the unalterability of the main proper¬ 
ties of the mock galaxies. 


5 PHOTOMETRIC REDSHIFT 
PERFORMANCE 

In this section, we investigate the behavior of the photomet¬ 
ric redshifts obtained from the mock catalogues described in 
section 3 for each next generation survey in terms of photo¬ 
metric redshift resolution and catastrophic outliers rate. 

5.1 Photometric Redshift Resolution 

In order to obtain a robust measurement of the mean pho¬ 
tometric redshift resolution, we have used the NMAD es- 
timato r following a similar approach as in iBrammer et al.l 
l|2008l) : iMolino et al.l ||2014| ). This estimator is known to be 
less sensitive to the behavior of the distribution. As such, 
the NMAD estimator is defined as: 

*0 >• /|A« - median{\z)\\ 

a nmad = 1-48 x median l -—-- 1 (1) 

where A z = Zb — z s 
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Figure 13. Stellar mass function for quiescent galaxies for six 
different redshift bins: 0.2 < z < 0.5 (top left panel), 0.5 < z < 
1.0 (top right panel), 1.0 < 2 < 1.5 (middle left panel), 1.5 < 2 < 
2.0 (middle right panel), 2.0 < 2 < 2.5 (bottom left panel) and 
2.5 < 2 < 3.0 (bottom right panel). The red diamonds refer to 
stellar mass counts of those galaxies classified as late-type (tj, > 3) 
in the original mock catalogue fit. The blue triangles refer to the 
stellar mass counts of those galaxies classified as late-type (tj, > 3) 
in the PhotReal post-processed mock catalogue for the Euclid- 
opt survey. The solid line and the gray area are the observational 
quiescent SM function obtained by Muzzin et al. 2013. 

In Fig 1171 we show the density maps of the photomet¬ 
ric redshift versus spectroscopic redshift for the three cases 
considered. We immediately notice several differences be¬ 
tween the considered surveys. In the case of the LSST, the 
lack of infrared data gets compensated with the depth of the 
optical data, being able to reach higher redshift and mag¬ 
nitude limits with a photometric redshift precision of 0.045 
for all the redshift range. We see some features at redshift 
~ 1.5, corresponding to the limit at which the spectrum of 
the objects is shifted to the infrared. The two cases con¬ 
sidered for the Euclid surveys sample well the optical and 
infrared spectrum providing a less biased and lower rates of 
catastrophic outliers than the other surveys. The expected 
accuracy for the Euclid pessimistic and optimistic case is 
cnmad = 0.070 and 0.0298 respectively. It is noticeable the 
higher dispersion found for the Euclid-Pes at lower redshift 
<0.4, which is prob ably due to the lack of UV data (i.e. 
ISchmidt fe Thormanll2013l ). 

We have now investigated the photometric redshift per¬ 
formance as a function of the magnitude. The results are 
collected in Table [l] and in Fig. 1181 we show the behavior 
of the photometric redshift bias as a function of magnitude. 
The parameter 77 refers to the fraction of outliers and it is 
defined in 1)5.21 Clearly, the number of bands covering a wide 
range of the spectrum has a direct impact in the photomet¬ 
ric redshift accuracy and photometric redshift bias. While 
the Euclid-Opt has photometric redshift precision of < 0.02 
down to H = 24 mag AB, we have an increase factor of 2-3 


Figure 14. Angular correlation function for the original mock 
catalogue (solid line) and the PhotReal post-processed mock cat¬ 
alogue (dotted line) for the sample limited to redshift 2 < 0.2. The 
shaded areas refer to 3<r errors. The color symbols and the dashed 
lines show the observational SDSS angular correlation function 
calculated by Wang et al. 2013. 
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Figure 15. Angular correlation function for the original mock 
catalogue (solid line) and the PhotReal post-processed mock cat¬ 
alogue (red dotted line) for the quiescent galaxies for six different 
redshift bins: 0.2 < 2 < 0.5 (top left panel), 0.5 < 2 < 1.0 (top 
right panel), 1.0 < 2 < 1.5 (middle left panel), 1.5 < 2 < 2.0 
(middle right panel), 2.0 < 2 < 2.5 (bottom left panel) and 
2.5 < 2 < 3.0 (bottom right panel). See the text for a defini¬ 
tion of quiescent galaxy for each sample. 
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Figure 16. Angular correlation function for the original mock 
catalogue (solid line) and the PhotReal post-processed mock cat¬ 
alogue (blue dotted line) for the star-forming galaxies for six dif¬ 
ferent redshift bins: 0.2 < z < 0.5 (top left panel), 0.5 < z < 1.0 
(top right panel), 1.0 < z < 1.5 (middle left panel), 1.5 < z < 2.0 
(middle right panel), 2.0 < z < 2.5 (bottom left panel) and 
2.5 < z < 3.0 (bottom right panel). See the text for a defini¬ 
tion of star-forming galaxy for each sample. 


for the Euclid-Pes and the LSST cases for the same range of 
magnitudes. The same effect is observed in the photometric 
redshift bias. 

We have also looked into the performance of the photo¬ 
metric redshifts for each different survey as a function of the 
redshift. The results are summarized in Tableland are dis¬ 
played in Fig. 1191 We find that the photometric redshift dis¬ 
persion is always at least a factor of 1.5 higher for the Euclid- 
Pessimistic than for the LSST. This result remarks that the 
presence of the infrared bands do not improve the photomet¬ 
ric redshift resolution at z < 1. By contrast, the improve¬ 
ment in the photometric redshift performance found in us¬ 
ing a deeper optical counterpart becomes remarkable when 
comparing the results of the two Euclid cases. The photo¬ 
metric redshift resolution of the Euclid-Optimistic is a factor 
of 2 to 5 better than in the Pessimistic case both in terms 
of photometric redshift accuracy and bias. In general, we 
find that the photometric redshift dispersion depends more 
strongly on m agnitude than on redshi ft , in agreement with 
other works (jHildebrandt et ah 20121: jPahlen et alj 120131 : 

dJl2C 


iRaichoor et al.ll2014l: iMolino et al. K 


5.2 Fraction of Outliers 

The term ’outliers’ refers to those galaxies that, due to the 
well-known color-redshift degeneracies, gets assigned a very 
similar color to a theoretical galaxy at a very different red¬ 
shift. Even if the Bayesian probability alleviates t his prob- 
lem due to the introduction of a Bayesian prior ( iBenit ez 


Figure IT. Density maps of spectroscopic redshift versus photo¬ 
metric redshift for each different survey considered for all the sim¬ 
ulation (LSST, upper right; Euclid pessimistic, upper left and Eu¬ 
clid optimistic, bottom left). For a definition of Euclid pessimistic 
and optimistic case see §3.3. The color scale is logarithmic and 
redder colors represents higher densities while bluer colors refer 
to less dens regions. The overall photometrical redshift precision 
for each survey is quoted in each panel. Individually, LSST quotes 
the highest photometric redshift dispersion in the whole redshift 
range due to the drop of the flux in the optical bands at high 
redshift. The Euclid pessimistic and optimistic cases sample well 
all the redshift range up to redshift 3, having the later an improve 
in overall photometric redshift precision of a factor of two with 
respect to the former. 


2000; Bellag amba et ahll2012l ), the existence of such galax¬ 
ies still remains and we need to quantify it. 

One necessary requirement for good quality photomet¬ 
ric redshifts is to have a small percentage of the total sample 
of ga laxies becoming ’outliers’ . We define fraction of outliers 
as in [Hildebrandt et all lj20ld ): 


|As| 

1 + z s 


> 0.15 


( 2 ) 


The results obtained for different surveys as a function 
of magnitude and redshift are collected in the third column 
of each survey in Tables [l] and [2] respectively. In Fig 1201 we 
show the outlier rate as a function of redshift for the different 
surveys. 

While the Euclid-Opt survey has a very lower rate of 
catastrophic outliers as a function of magnitude, being less 
than 1% down to H = 23 mag, the LSST and Euclid-Pes 
shows a higher rate of outliers between 3-12%. This sepa¬ 
ration illustrates that the multiple (> 12) band surveys are 
much more effective in decreasing the rate of catastrophic 
outliers than surveys with fewer bands. 

Looking at the same rate as a function of redshift, we 
see a milder dependence of the outlier rate and redshift. The 
rate of outliers remains close to ~l-2% for the Euc-Opt, ~ 
2-10% for the LSST and~ 10-20% for the Euc-Pes survey. 
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Table 1. Photometric redshift comparison for different surveys as a a function of magnitude 





LSST 



Euclid — P 



Euclid — O 


mag 

Az 

ONMAD 

V 

Az 

°NMAD 

V 

Az 

°NMAD 

V 

16.0 - 

16.5 

0.0276 

0.0333 

15.3281 

0.0239 

0.0342 

19.2475 

0.0040 

0.0137 

0.2043 

16.5 - 

17.0 

0.0242 

0.0312 

13.2013 

0.0206 

0.0324 

15.9318 

0.0041 

0.0137 

0.1542 

17.0 - 

17.5 

0.0199 

0.0290 

10.9262 

0.0157 

0.0300 

11.9770 

0.0030 

0.0143 

0.1742 

17.5 - 

18.0 

0.0176 

0.0287 

9.4708 

0.0113 

0.0309 

9.8102 

0.0033 

0.0152 

0.1491 

18.0 - 

18.5 

0.0156 

0.0277 

7.7501 

0.0063 

0.0318 

8.1134 

0.0023 

0.0152 

0.1668 

18.5 - 

19.0 

0.0124 

0.0273 

6.6554 

0.0016 

0.0322 

7.2817 

0.0014 

0.0153 

0.1645 

19.0 - 

19.5 

0.0096 

0.0268 

5.5151 

-0.0015 

0.0315 

7.1102 

0.0008 

0.0151 

0.1601 

19.5 - 

20.0 

0.0067 

0.0260 

4.6757 

-0.0031 

0.0300 

7.2907 

0.0007 

0.0150 

0.1496 

20.0 - 

20.5 

0.0046 

0.0251 

3.8189 

-0.0034 

0.0287 

7.4268 

0.0009 

0.0150 

0.1319 

20.5 - 

21.0 

0.0029 

0.0242 

3.1238 

-0.0037 

0.0277 

7.5001 

0.0011 

0.0150 

0.1296 

21.0 - 

21.5 

0.0015 

0.0237 

2.6430 

-0.0040 

0.0278 

7.4401 

0.0011 

0.0150 

0.1304 

21.5 - 

22.0 

0.0006 

0.0235 

2.2978 

-0.0046 

0.0291 

7.2443 

0.0010 

0.0154 

0.1564 

22.0 - 

22.5 

-0.0001 

0.0236 

2.1691 

-0.0061 

0.0323 

7.1957 

0.0006 

0.0163 

0.2341 

22.5 - 

23.0 

-0.0005 

0.0238 

2.3066 

-0.0089 

0.0396 

8.9712 

-0.0001 

0.0182 

0.4693 

23.0 - 

23.5 

-0.0008 

0.0242 

2.7025 

-0.0123 

0.0546 

15.2487 

-0.0006 

0.0216 

1.2641 

23.5 - 

24.0 

-0.0007 

0.0252 

3.5816 

-0.0208 

0.0791 

26.1357 

-0.0006 

0.0277 

3.5440 

24.0 - 

24.5 

-0.0007 

0.0278 

5.2341 







24.5 - 

25.0 

-0.0005 

0.0327 

7.5035 







25.0 - 

25.5 

-0.0051 

0.0441 

12.9205 







25.5 - 

26.0 

0.0031 

0.0657 

23.2505 








Table 2. Photometric redshift comparison for different surveys as a a function of redshift 





LSST 



Euclid — P 



Euclid — O 


z 

Az 

CNMAD 

V 

Az 

^NMAD 

V 

Az 

ff NMAD 

V 

0.0- 

0.5 

0.0057 

0.0302 

5.8276 

o.oou 

0.0457 

13.6601 

0.0009 

0.0214 

1.4644 

0.5 - 

1.0 

-0.0021 

0.0237 

2.0035 

-0.0225 

0.0449 

17.8672 

-0.0007 

0.0174 

1.0415 

1.0 - 

1.5 

0.0143 

0.0300 

10.3601 

0.0023 

0.0414 

12.4176 

0.0013 

0.0171 

0.8455 

1.5 - 

2.0 

-0.0199 

0.0452 

11.0907 

0.0089 

0.0564 

15.1593 

-0.0002 

0.0291 

2.4469 

2.0- 

2.5 

-0.0409 

0.0401 

12.4108 

-0.0429 

0.0581 

14.7191 

-0.0017 

0.0285 

3.0074 

2.5 - 

3.0 

-0.0181 

0.0229 

7.2210 

-0.1783 

0.0577 

20.8091 

-0.0048 

0.0166 

1.8212 


As before, the presence of the deep infrared data and the 
combination with deep optical data is a strong combination 
to decrease the rate of catastrophic outliers. 


6 PHOTOMETRIC REDSHIFT QUALITY 

CUTS 

As we have illustrated in the previous section, the existence 
of outliers in a photometric survey is unavoidable but can be 
quantified and keep under a moderate threshold. In addition, 
we can use one main advantage of some photometric red¬ 
shift packages such as BPZ2.0, which is the fact that all the 
galaxies have assigned a redshift probability function. The 
study of this function provides a measure of the quality of 
the photometric redshift: very well behaved photometry red- 
shifts will probably have a well defined narrow peak around 
the spectroscopic redshift of the galaxy, whereas galaxies 
with degeneracies will be assigned a multi peak probability 
dist ributi on. 

iBem'ted J2000l l introduced the odds parameter: a mea¬ 
surement of the quality of a photometric redshift. This 
odds parameter is obtained by integrating the redshift 
probability function of the galaxy within 2<r of its more 
probable value. These parameters has been widely used 
to select the best quality photometric redshi ft sam ple to 
study a number of different features (e.g. ICoe et ahl l200d : 


Benitez et al. 2( 

)09a,b; Bellagainba et al. 2012|; Molino et al. 

20141: Ascaso et 

al. 2015(') and allows to improve the overall 


photometric redshift precision of the sample by previously 
modeling the selection bias. 


In this section, we discuss the benefit obtained in terms 
of the improvement of the photometric redshift accuracy and 
the decrement of the outlier rates when using this parameter 
to select the best quality photometric redshift sample for a 
particular survey. 


6.1 Improvement in the photometric redshift 
resolution 

In order to illustrate the increase in photometric redshift 
resolution by effectuating a cut in the odds parameter, we 
plot in Fig. 1211 the photometric redshift accuracy as a func¬ 
tion of the odds parameter for each of the three different sur¬ 
veys considered in this study. As expected, we see a decrease 
of the photo-z dispersion as the odds increase for all three 
surveys. Indeed, we can select a subsample of the accuracy 
desired by performing an accurate odds cut. For instance, 
while the overall accuracy of LSST is 0.045, we could obtain 
a subsample of twice this accuracy (i.e. ijnmad ~ 0.0279) 
by selecting the 25% highest odds subsample. 

We also illustrate the effect of this selection by plotting 
the analogous plot to the Fig. 1171 in Fig. 1221 the behavior of 
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Figure 18. Density maps of magnitude versus the photometric 
bias for the three different survey considered in this work (Euclid 
pessimistic, top left; Euclid optimistic, bottom left and LSST, top 
right).The color scale are the same as in Fig |17l The axis scale are 
the same for all the panels. The behavior of the photometric red- 
shifts for the LSST and both Euclid cases show larger dispersions 
and biases down to deeper magnitudes, having the Euclid Opti¬ 
mistic case a significantly smaller dispersion than the Pessimistic 
case and the LSST. 


the photometric redshift as a function of spectroscopic red- 
shift for the 50% best quality photometric redshifts for each 
of the surveys considered. The number quoted in the legend 
of each of the panels refers to the overall precision achieved, 
being this a factor between 1.5 and 3 with respect to the 
overall sample. In particular, a large fraction of outliers, de¬ 
parting from the one-to-one line in Fig. 1171 disappears by 
considering only galaxies with high odds. 


6.2 Decrease of outlier rates 

As mentioned in the previous section, in Fig. 1201 we have il¬ 
lustrated the relation between the outlier rates as a function 
of redshift with solid lines for each of the different surveys. 
In this Figure, we have also plotted the outlier rates result¬ 
ing for the 50% best photometric redshift quality subsample 
for every survey. We plainly observe a dramatic descend of 
the outlier rate for each survey. The new rates are decreased 
by a factor of 3 to 10, putting in evidence the virtues of this 
selection. 

When working in a particular science topic, one needs 
to understand the nature of any possible selection based on 
the odds parameter and quantify the possible bias intro¬ 
duced on the final scientific results. Once this effect is taken 
into account, the improvement in photometric precision and 
the decrease of outliers rate can be fully exploited and we 
encourage researchers to exploit this feature. 



0 12 3 


Figure 19. Density maps of redshift versus the photometric bias 
for the LSST, upper right; Euclid pessimistic, upper left and Eu¬ 
clid optimistic, bottom left. The color scale are the same as in Fig 
1171 The axis scale are the same for all the panels. The behavior for 
the LSST and Euclid-Pessimistic case becomes similar (less than 
a factor of 2) up to redshift 1.5. At higher redshift, the bias for the 
Euclid-Pessimistic case, than the Euclid Optimistic case, which 
obtains an overall good performance (ctnmad ~ 0.02 — 0.4) and 
negligible photometric redshift bias up to redshift 3. 

7 DISCUSSION AND CONCLUSIONS 

The advent of several large deep surveys is imminent and 
the need of precise tools to describe accurately their proper¬ 
ties is pressing. In this work, we have developed mock cat¬ 
alogues mimicking realistically two of the up-coming gen¬ 
eration stage IV surveys: LSST and Euclid (considering 
two cases for the optical counterpart, optimistic and pes¬ 
simistic). We have confirmed that the main properties of 
this new mock reflect reality in terms of photometry, photo¬ 
metric errors, stellar masses, spectral types and photometric 
redshifts by comparing with observable relations. We have 
also provided predictions in terms of photometric redshift 
performance and selection cuts from the same basis for each 
survey. The mocks will become publicly availabk{f] and they 
will become extremely useful for making internal compar¬ 
isons between the different surveys. In the appendix, we in¬ 
clude a description of the content of the mock catalogues. 

Each of these surveys will provide output data observed 
with very different strategies. In this work, we have disen¬ 
tangled some of the advantages of using each of these dif¬ 
ferent datasets in terms of photometric redshift resolution. 
First, the number of bands used to sample the optical and/or 
infrared spectrum is directly related to the level of photo¬ 
metric redshift accuracy and photometric redshift bias as 
confirmed by other works (e.g. iBenftez et al.ll2009al and ref¬ 
erences herein). On the other hand, deep broad-band opti- 


http: //photmocks. obspm.fr/ 
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Figure 20. Fraction of outliers as a function of the redshift for 
the three different survey considered in this work (LSST, red line; 
Euclid pessimistic, blue line and Euclid optimistic, green line). 
The solid lines refer to the overall population for each survey 
while the dotted lines refer to the 50% of the population with 
highest odds for each different survey. In general, the fraction of 
outliers increases gently with redshift for all the surveys consid¬ 
ered here, obtaining the lower rates for Euclid Optimistic, Euclid 
Pessimistic, and LSST in increasing order. If we consider the best 
50% quality odds sample, the outliers rate decreases more than a 
factor of 10. 

cal surveys allow us to sample the luminosity function in 
the optical range even if with worse photometric redshift 
resolution. They are also an excellent complement for deep 
infrared surveys, particularly to dramatically decrease the 
rate of outliers and improve the photometric redshift accu¬ 
racy. Finally, the deep infrared surveys provide an excellent 
photometric redshift accuracy in a wide range of redshift 
(z < 3) with a moderate rate of catastrophic outliers. Thus, 
the future combination of all these strategies can lead to 
large amount of ancillary data which will provide excellent 
resolution for many scientific purposes. These mocks will 
become extremely useful for making internal comparisons 
between different surveys. 

As stated in the introduction, a separate work (Ascaso 
et al. in prep) will be fully devoted to explore the photo¬ 
metric redshift expectations for J-PAS, a multiple narrow- 
band survey. We have then examined the performance of 
the photometric redshifts for a new ’super-survey’ consist¬ 
ing in combining all the expected data for these surveys 
once they are made public and we are allowed to combine 
it. This ’super survey’ will consist then in the combination of 
Euclid+LSST+DES+J-PAS and the results should account 
for all the advantages described through this work. In Fig |23l 
we plot the density map of the photometric redshift versus 
the spectroscopic redshift for this ’super survey’. Let’s note 
that these three surveys will not overlap in all parts of the 
sky and their overlap will also be different in the southern 
and northern hemispheres. 



Odds 

Figure 21. Odds photometric redshift quality parameter versus 
photometric redshift dispersion for the three different survey con¬ 
sidered in this work (LSST, red dotted line; Euclid pessimistic, 
blue dashed line and Euclid optimistic, green dotted-dashed line). 
Performing a safe quality odds parameter in each of these surveys 
can decrease the photometric redshift dispersion significantly. 



Z s 


Figure 22. Density maps of spectroscopic redshift versus photo¬ 
metric redshift for the 50% best quality photometric redshift for 
each different survey considered for all the simulation (LSST, up¬ 
per right; Euclid pessimistic, bottom left and Euclid optimistic, 
bottom right). The color scale are the same as in Fig | 1TI The over¬ 
all photometrical redshift precision for each survey is quoted in 
each panel. The photometric redshift resolution increases between 
a factor of 1.5 and 3 for all surveys by selecting a best quality 
photometric redshift subsample, being particularly affected the 
density regions of outliers. 
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Figure 23. Density map of spectroscopic redshift ver¬ 
sus photometric redshift for a future ’super survey’ being 
Euclid+LSST-|-DES-|-J-PAS. The color scale are the same as in 
Fig |17l The behavior of the photometric redshift when the ancil¬ 
lary data of some of the most important next generation survey 
is available would be excellent. The catastrophic outlier rate will 
decrease dramatically and the 0 < z < 3 redshift range will be 
measured with an accuracy of 0.0095 without imposing any odds 
quality cut. 




H mag [AB] z 

Figure 24. Catastrophic outliers rate as function of H-band (left 
plot) and redshift (right plot) for a future ’super survey’ being 
Euclid+LSST-|-DES-|-J-PAS. The catastrophic outlier rate will 
never be higher than 3% in all magnitude and redshift ranges. 


The expected overall photometric redshift resolution for 
all the redshift range within 0 < z < 3 and i < 27.5 AB will 
be of cy nmad = 0.0135, with almost negligible photometric 
redshift bias. As for the rate of catastrophic outliers, we plot 
in Fig. 1241 the expected outlier rate for this ’super-survey’ as 
a function of magnitude and redshift. As we see, this survey 
will keep the outliers rate below 1% up to redshift ~2 and 
below 3% for all the redshift and magnitude range. 

Hence, this future ’super survey’ will support almost all 
the scientific cases for the community, becoming a superior 


machine comparable to a full spectroscopic survey ranging 
a very large range of redshift and magnitude. 

In a forthcoming paper (Ascaso et al., in prep), we will 
continuos this work within the frame of Apples to Apples, 
A 2 to explore and compare the expected selection function 
for cluster and group detections for the three surveys con¬ 
sidered here. We will also publish a third paper providing 
cosmological constraints from these cluster counts for dif¬ 
ferent surveys. These series of works will help setting light 
on the expected predictions of these surveys and the mocks 
will be useful to set requirements to the science derived from 
them. 
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APPENDIX A: MOCK CATALOGUES 
DESCRIPTION. 

In this appendix we include the description of the content 
of the content of the galaxy mock catalogues in detail. The 
three mock catalogues are 500 deg 2 wide and the magnitudes 
limits for the different bands are: 

• Euclid: magnitude limit of H=24, J=24 and Y=24, 
down to 5a, (ILaureiis et al]|2011f ). 

• LSST: magnitude limit of u=26.1 , g=27.4, r =27.5, 
i=26.8, z=26.1 and y=24.9, down to 5a, (llvezic et al.ll2008l) . 

• DES: magnitude limit of g=25.2, r=24.8, i=24.0, 
z=23.4 and y=21.7, down to lOcr, dMohr et alj|20l3 l. 
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Table Al. Content of the Euclid Optimistic Mock Catalogue 


PARAMETER 

DESCRIPTION 

DHaloID* 

ID of the dark matter halo in the Millennium Simulation 

Galaxy ID* 

ID of the galaxy 

is.central* 

= 1 if the galaxy is central, 0= satellite 

MHalo* 

Mass of the host halo (Mq) 

MStars.tot* 

Total Stellar mass of the galaxy (Mq) 

RA* 

Right Ascension in decimal degrees [J2000] 

DEC* 

Declination in decimal degrees [J2000] 

zs* 

Spectroscopic redshift 

SDSS_u_or* 

SDSS.u original magnitude 

SDSS_g_or* 

SDSS_g original magnitude 

SDSS_r_or* 

SDSS_r original magnitude 

SDSS_i_or* 

SDSS_i original magnitude 

SDSS.z.or* 

SDSS_z original magnitude 

DES_g_or* 

DES_g original magnitude 

DES_r_or* 

DES_r original magnitude 

DES.i.or* 

DES_i original magnitude 

DES_z_or* 

DES_z original magnitude 

DES_y_or* 

DES_y original magnitude 

Euclid_H_or* 

Euclid_H original magnitude 

Euclid_J_or* 

Euclid.J original magnitude 

Euclid_Y_or* 

Euclid_Y original magnitude 

zb 

BPZ2.0 most likely redshift 

zb_min 

Lower limit (95p confidence) 

zb_max 

Upper limit (95p confidence) 

tb 

BPZ2.0 most likely spectral type 

Odds 

P(z) integrated within zb - 2 &NMAD and zb + 2 cr^MAD 

Chi2 

Poorness of BPZ fit: observed vs. model fluxes 

StelLMass 

PhotRealStellar Mass (loglO(MQ)) 

DES_g, e_DES_g 

DES_g PhotReal magnitude uncertainty 

DES_r, e_DES_r 

DES_r PhotReal magnitude & uncertainty 

DES J, e.DESJ 

DES_i PhotReal magnitude & uncertainty 

DESj, e_DES_z 

DES_z PhotReal magnitude & uncertainty 

DES.y, e_DES_y 

DES_y PhotReal magnitude & uncertainty 

LSST.u, e_LSST_u 

LSST.u PhotReal magnitude & uncertainty 

LSST_g, e_LSST_g 

LSST_g PhotReal magnitude & uncertainty 

LSST_r, e_LSST_r 

LSST_r PhotReal magnitude & uncertainty 

LSSTJ, e.LSSTj 

LSSTJ PhotReal magnitude & uncertainty 

LSSTj, e_LSSTj 

LSST.Z PhotReal magnitude & uncertainty 

LSST.y, e_LSST_y 

LSST.y PhotReal magnitude & uncertainty 

Euclid_H, e_Euclid_H 

Euclid_H PhotReal magnitude & uncertainty 

Euclid.J, e_Euclid_J 

Euclid.J PhotReal magnitude & uncertainty 

Euclid.Y, e_Euclid_Y 

Euclid.Y PhotReal magnitude & uncertainty 


as in the original Merson et al. mock catalogue 




I. Apples to apples A 2 : photo-z predictions for next-generation surveys 
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Table A2. Content of the Euclid Pessimistic Mock Catalogue 


PARAMETER 

DESCRIPTION 

DHaloID* 

ID of the dark matter halo in the Millennium Simulation 

Galaxy ID* 

ID of the galaxy 

is.central* 

= 1 if the galaxy is central, 0= satellite 

MHalo* 

Mass of the host halo (Mq) 

MStars.tot* 

Total Stellar mass of the galaxy (Mq) 

RA* 

Right Ascension in decimal degrees [J2000] 

DEC* 

Declination in decimal degrees [J2000] 

zs* 

Spectroscopic redshift 

SDSS_u_or* 

SDSS.u original magnitude 

SDSS-g-or* 

SDSS_g original magnitude 

SDSS_r_or* 

SDSS_r original magnitude 

SDSS_i_or* 

SDSS_i original magnitude 

SDSS.z.or* 

SDSS_z original magnitude 

DES_g_or* 

DES_g original magnitude 

DES-r-or* 

DES_r original magnitude 

DES.i.or* 

DES_i original magnitude 

DES_z_or* 

DES_z original magnitude 

DES_y_or* 

DES.y original magnitude 

Euclid_H_or* 

Euclid_H original magnitude 

Euclid_J_or* 

Euclid.J original magnitude 

Euclid_Y_or* 

Euclid.Y original magnitude 

zb 

BPZ2.0 most likely redshift 

zb_min 

Lower limit (95p confidence) 

zb_max 

Upper limit (95p confidence) 

tb 

BPZ2.0 most likely spectral type 

Odds 

P(z) integrated within zb - 2 &NMAD and zb + 2 cr^MAD 

Chi2 

Poorness of BPZ fit: observed vs. model fluxes 

StelLMass 

PhotReal Stellar Mass (loglO(MQ)) 

DES_g, e_DES_g 

DES_g PhotReal magnitude & uncertainty 

DES_r, e_DES_r 

DES_r PhotReal magnitude & uncertainty 

DES J, e.DESJ 

DES_i PhotReal magnitude &; uncertainty 

DESj, e_DES_z 

DES_z PhotReal magnitude &; uncertainty 

DES.y, e_DES_y 

DES.y PhotReal magnitude & uncertainty 

Euclid_H, e_Euclid_H 

Euclid_H PhotReal magnitude uncertainty 

Euclid.J, e_Euclid_J 

Euclid.J PhotReal magnitude & uncertainty 

Euclid.Y, e_Euclid_Y 

Euclid.Y PhotReal magnitude & uncertainty 


as in the original Merson et al. mock catalogue 
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Table A3. Content of the LSST Mock Catalogue 


PARAMETER 

DESCRIPTION 

DHaloID* 

ID of the dark matter halo in the Millennium Simulation 

Galaxy ID* 

ID of the galaxy 

is.central* 

=1 if the galaxy is central, 0= satellite 

MHalo* 

Mass of the host halo (Mq) 

MStars.tot* 

Total Stellar mass of the galaxy (Mq) 

RA* 

Right Ascension in decimal degrees [J2000] 

DEC* 

Declination in decimal degrees [J2000] 

zs* 

Spectroscopic redshift 

SDSS_u_or* 

SDSS.u original magnitude 

SDSS-g-or* 

SDSS_g original magnitude 

SDSS_r_or* 

SDSS_r original magnitude 

SDSS_i_or* 

SDSS_i original magnitude 

SDSS-Z.or* 

SDSS_z original magnitude 

DES-g-or* 

DES^g original magnitude 

DES_r_or* 

DES_r original magnitude 

DES-Lor* 

DES J original magnitude 

DES_z_or* 

DES_z original magnitude 

DES_y_or* 

DES_y original magnitude 

zb 

BPZ2.0 most likely redshift 

zb_min 

Lower limit (95p confidence) 

zb_max 

Upper limit (95p confidence) 

tb 

BPZ2.0 most likely spectral type 

Odds 

P(z) integrated within zb - 2a n mad and zb + 2a nm AD 

Chi2 

Poorness of BPZ fit: observed vs. model fluxes 

StelLMass 

PhotReal Stellar Mass (loglO(M q)) 

LSST.u, e_LSST_u 

LSST.u PhotReal magnitude & uncertainty 

LSST_g, e_LSST_g 

LSST_g PhotReal magnitude &; uncertainty 

LSST_r, e_LSST_r 

LSST_r PhotReal magnitude & uncertainty 

LSSTJ, e_LSST j 

LSSTJ PhotReal magnitude & uncertainty 

LSST_z, e_LSST_z 

LSST.Z PhotReal magnitude & uncertainty 

LSST.y, e_LSST_y 

LSST.y PhotReal magnitude & uncertainty 


as in the original Merson et al. mock catalogue 




